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We study one-dimensional sideband cooling of Cesium atoms strongly confined in a far-detuned 
optical lattice. The Lamb-Dicke regime is achieved in the lattice direction whereas the transverse 
confinement is much weaker. The employed sideband cooling method, first studied by Vuletic et 
al.Q, uses Raman transitions between Zeeman levels and produces a spin-polarized sample. We 
present a detailed study of this cooling method and investigate the role of elastic collisions in the 
system. We accumulate 83(5)% of the atoms in the vibrational ground state of the strongly confined 
motion, and elastic collisions cool the transverse motion to a temperature of 2.8 fiK—0.7 hui osc /kB, 
where uj oec is the oscillation frequency in the strongly confined direction. The sample then approaches 
the regime of a quasi-2D cold gas. We analyze the limits of this cooling method and propose a 
dynamical change of the trapping potential as a mean of cooling the atomic sample to still lower 
temperatures. Measurements of the rate of thermalization between the weakly and strongly confined 
degrees of freedom are compatible with the zero energy scattering resonance observed previously in 
| weak 3D traps. For the explored temperature range the measurements agree with recent calculations 

^•f^ of quasi-2D collisions Q. Transparent analytical models reproduce the expected behavior for feT 2> 

■ Juj osc and also for k^T <C TwJosc where the 2D features are prominent. 

Q , PACS numbers: 32.80Pj,34.50-s 



INTRODUCTION 



The physics of systems with reduced spatial dimensionality has attracted a great deal of interest. 
^ \ These systems, where one (two) degree(s) of freedom are confined to the quantum ground state of 
H ' motion, have properties which can markedly differ from those of 3D systems. The case of a 2D 
^ | Bose gas has been explored theoretically and experimentally with atomic hydrogen adsorbed on liquid 
• <-h . helium (see || |J] ) • Advances on laser cooling and Bose-Einstein condensation of atomic gasesJH] 
have opened an interesting possibility to create new ID and 2D quantum degenerate systems W. 
In a recent paper Petrov et al. showed that collisional properties of cold gases can be drastically 
modified by strongly confining the motion of particles in one direction^ : for instance, in a gas with 
negative scattering length a, the mean-field interaction can switch sign under variations of the strong 
confinement. This is of particular interest for the case of cesium. Due to a very large and negative 
scattering length a = — 138 nm in the \F = 3, to — 3) state a cesium condensate in a weakly 
confining trap would be unstable for a few tens of atoms |J. The characteristic size of the quantum 
ground state in the strongly confined direction is another important length in this problem. For a 
harmonic potential with oscillation frequency uj osc this length is la = y/Ti/2mu) osc . Ref. || predicts 
that for a sufficiently large value of the ratio a/lo the mean- field interaction becomes positive allowing 
the formation of a stable condensate. A second interest in a strong confinement in one direction is the 
possibility to use a very efficient optical cooling method, the sideband cooling [h], [nj |l2|, [u] . 
This method requires Iq < A, (Lamb-Dicke regime) where A is the optical wavelength of the cooling 
transition. For instance in ref.|T^], 92(5)% of the atoms have been cooled to the ground state in the 
tightly confined direction. Several trap configurations for achieving 2D gases are under investigation. 
Many of them are based on far-off resonance optical dipole traps pa] and evanescent fields near the 
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surface of a dielectric [17]] . Another common approach is to confine atoms in micro-wells of a ID 
far-detuned optical lattice [l2|, [l^] . 

In this article, we present an implementation of sideband cooling of cesium atoms in an intensity 
lattice created at the intersection of two far detuned YAG laser beams. The cooling method, first 
studied in rcf.Jjj], uses Raman transitions between Zeeman substates of \F = 3) and produces an atomic 
sample polarized in \F — 3, to — 3). We give a detailed account of cooling in our geometry which is 
different from that in ref.jjj. We recall the cooling principle, analyze the role of the beam polarizations 
and discuss the steady-state of the one-dimensional cooling. The intensity lattice is vertical but atoms 
are also weakly confined in the horizontal plane owing to the YAG gaussian beam intensity profiles. 
Elastic collisions couple the vertical and horizontal degrees of freedom so that ID sideband cooling 
efficiently cools in the three directions [Q . A horizontal temperature — 2.8 fiK— 0.7hui osc /kB has 
been obtained, which already corresponds to the condition of a quasi two-dimensional gas. We discuss 
the limits of this cooling. In particular, we show that the cooling efficiency becomes exponentially 
small for k^T < huj osc and propose a dynamical method to reach a lower ratio k^T /Tiuj OS c- Finally 
we address the question of the influence of the strong confinement on collisional properties. We 
measure the rates of thermalization between the strongly and weakly confined degrees of freedom. 
In the explored temperature range from 4 to 20 /iK, our results are compatible with the zero energy 
scattering resonance previously observed in weak 3D traps jl9], . The experimental data also agree 
with calculations of ref . Q which take into account the quantum character of the particle motion in 
the strongly confined direction. We propose simple analytical models that reproduce the expected 
behaviors for k^T ^> hu> osc , and also for k^T <C hu osc where the 2D character of the particle motion 
is prominent. 



II. THE FAR DETUNED TRAP 

Our ID intensity lattice is similar to the one described in fl2f . It is produced at the crossing of two 
beams of a YAG laser (1.06 /im), propagating in a vertical plane and making an angle 6*yag = 52° 
with the horizontal direction. Both beams have linear polarisations as presented in fig. [j]. Each beam 
has a waist of 100 /im and a power of about 5 W. The YAG laser is far detuned to the red of the 
Dl and D2 transitions of cesium at 852 nm and 894 nm. The atoms experience a dipole potential 
proportional to the laser intensity that confines the atoms in regions of maximum intensity. For 
horizontal polarizations, the interference of the two beams creates a vertical intensity lattice with a 
period of 665 nm, and the dipole potential is modulated in the vertical direction. The horizontal 
confinement is provided by the gaussian shape of the YAG beams. The total trap depth is about 
140 /iK. Tunneling between different micro-wells is totally negligible for temperatures smaller than 
50 /iK and the atoms are confined in independent micro-traps. In the central micro-trap, the vertical 
oscillation frequency is lu osc /2tt = 80 kHz and the horizontal oscillation frequencies are u x /2n =175 Hz 
and ixJy/2-K = 140 Hz. The spontaneous emission rate at the bottom of the trap is about 3s _1 . It leads 
to a heating rate of 0.38 /iK/s that we can neglect in the following experiments. This trap is loaded 



from a magneto-optical trap as explained in 12 . About 2 x 10 5 atoms are trapped in a gaussian 
cloud of rms dimensions a x = 31 /im, a y = 39/im, a z = 60/iin, thus populating about 200 horizontal 
micro-traps. Their temperature is ~ 20 /iK, which for the vertical motion corresponds to a mean 
vibrational number (n) = 5.8. The vertical oscillation frequency depends on the position in the trap 
and varies by 15% within the spatial extension of the cloud. 

Using a charged-coupled-device camera, we perform two-dimensional absorption images of the 
atomic cloud. The probe beam is horizontal and makes an angle of 45° with the plane of the YAG 
beams. First, an image taken just after switching off the YAG beams gives a measurement of the 
trap size. Second, an image taken after a free expansion time of the cloud gives access to the velocity 
distribution in the vertical and horizontal directions in the plane of the camera. This time of flight 
method is detailed in EJ . 

The Lamb-Dicke parameter is rj = ^uj ycc /uj osc = 0.16 where w rcc = hk 2 /2m is the recoil frequency 
associated with the D2 transition. This low value enables us to cool the atoms to the ground state of 
motion in the vertical direction by using sideband cooling methods jl0| O, O, ^2) . 
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FIG. 1: (a) : Laser beam configuration used in sideband cooling of cesium. Atoms are confined in independent 
horizontal microwells resulting from the interference between the two YAG laser beams propagating in the yz 
plane. A weak magnetic field is applied along x. A repumping beam tuned to \Qs\/2,F = 3) — > \&Pa/2,F' = 2) 
propagates in the xy plane and makes an angle of 45° with the x axis. The polarization of one YAG beam is 
linear along x and the polarization of the other YAG beam (MB) makes a small angle a with x. (b) : Trapping 
potential along z. For more visibility, the period of the lattice (665nm) has been increased by a factor 40. 



III. COOLING METHOD 



A. Principle 



We describe here the main elements of the sideband cooling of the vertical motion as developed by 
Vuletic et al.jlj. This cooling uses only Zeeman sublevels of the Cs F = 3 hyperfine ground state and 
produces a sample polarized in \F = 3, m — 3). A magnetic field is applied along the x axis such that 
the Zeeman sublevels of F = 3 acquire different energy shifts. A repumping laser which is resonant 
with the transition |6s]/2, F = 3) — > |6p 3 / 2 , F' — 2) and which has a polarization with only it and a + 
components, provides a finite lincwidth of all Zeeman substates of F = 3 except for \F = 3, m = 3). 
Zeeman states with m differing by ±1 are coupled by a Raman process such that the Raman transition 
|m) — > \m — 1) is resonant with the transition that decreases the vibrational quantum number by one. 
If the Raman coupling and the intensity of the repumping laser are sufficiently small, the motional 
sidebands \m — 3, n) — * |m = 2, n') are well resolved and the state \m = 3,n = 0) is a quasi-dark 
state. Its lifetime is limited only by off-resonant Raman transfer and is much longer than that of 
the other states. In the Lamb-Dicke regime the vibrational energy of the atoms is much larger than 
the recoil energy. Therefore, the Raman transition to a Zeeman level with m < 3 followed by optical 
repumping to m = 3, leads on average to a decrease of the atoms' energy. Thus, atoms are cooled 
by repetition of such cycles and they accumulate in the quasi-dark state \m = 3,n = 0). A weak 
intensity laser tuned to the \F — 4) — > \F' — 4) transition prevents atoms from accumulating in the 
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FIG. 2: Number of atoms transferred to the hyperfine ground state \F — 4) by a two-photon Raman pulse as 
a function of the detuning of the two-photon transition from 9.192 631 770 GHz. The Raman beams propagate 
in the same direction in order to avoid change of vibrational state during the transfer. The labels indicate the 
various \F = 3, m) — » |_F = 4, m') transitions. Dashed line : unpolarized atoms. Solid line : atoms polarized 
by a 3.6ms repumping laser pulse with an intensity of 0.3/ sa t. The resulting polarisation exceeds 90%. 



\F — 4) state after spontaneous emission. The magnetic field is chosen so that the energy of |m, n) 
equals the energy of |m — 1, n — 1). The Raman transition is then resonant if the two Raman beams 
have the same frequency, and in contrast to our previous work|ll], |l^| we simply use the YAG beams 
themselves to induce the Raman transitions. In order to introduce coupling between different Zeeman 
levels, the p olariz ation of one of the beams is slightly modified from the linear horizontal polarization 
(see section |III Cf ). 



B. Repumping beam 



The repumping beam tuned to the |.F = 3) — > \F' = 2) transition of the D2 line propagates in 
the xy plane and makes an angle of 45° with the x axis. With respect to the axis of the magnetic 
field , the total intensity of the beam / splits in 1/3 for the it polarization (i.e. along the axis of the 
magnetic field), and 21/3 for the er + polarization. 

The population of the atoms in each of the Zeeman sublevels of F = 3 is measured via a frequency 
selective Raman transitions to \F = 4). For this purpose we use two Raman beams with frequencies 
differing by about 9.2 GHz as explained in p^ j. 

With this method we optimize the parameters of the repumping beam in the absence of Raman 
couplin g bet ween Zeeman sublevels i.e. with parallel linear polarizations of the YAG beams (see 
section [II C). We can measure both the polarisation time constant and the steady-state polarization. 
With our experimental precision, the equilibrium population in \m — 3) deduced from the data in fig. 
||is between 90% and 100%. We then deduce an intensity for the er_ component of the repumping laser 
smaller than 3% of the n intensity. The measured 1/e polarization time of an initially unpolarized 
sample is 10 fj,sxl sat /l at small intensities. The calculated lifetime T' -1 of the state |m = 2) is 
0.8 nslgat/l- The optical density of our sample is about 2 and we have seen no influence of the 
number of atoms on the polarization time and on the equilibrium state. This means that we see no 
effect of reabsorption of photons on the polarisation. 



C. Polarization of the YAG beams 



The YAG detuning is much larger than the hyperfine splittings of the 6si/2 and 6733/2 states. 
Hence, if the two YAG beams are linearly polarized along x, there is no coupling between different 
Zeeman sub-levels. The light-shift operator is scalar (proportional to the identity operator) and 
all Zeeman substates see the same YAG potential |l5|. Raman coupling between different Zeeman 
levels is obtained by introducing, for one of the trapping beams (called MB), a small component of 
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FIG. 3: Raman coupling V(z) ( dashed line) between \F = 3, m = 3) and \F = 3,m = 2) as a function of the 
vertical position of the atom. Figure (a) corresponds to a linear polarization of the MB making an angle of 20° 
with the horizontal direction and figure (b ) corresponds to a polarization of the MB with the same proportion of 
horizontal polarization but with a relative phase of n/2 between the two polarizations. The solid lines depicts 
the dipole potential seen by the atoms, the potential being chosen at the bottom of the micro-wells. For better 
visibility, we have plotted 80 x V(z). 



polarization X° orthogonal to the horizontal one. All Zeeman sublevels still share a common potential 
energy and levels |m,n) and \m — l,n — 1) are degenerate for any m £ [—2,3] and vibrational level 
n > 0. Raman transfers enable a change of the vertical motion. This corresponds to absorption 
from one beam and stimulated emission into the other beam. Due to an interference between the two 
different processes, at the bottom of a micro-trap the Raman coupling depends on the relative phase 
between the component of horizontal polarization and the component along X° . The coupling is even 
in z for a phase difference of ir/2 ( most elliptical polarisation of the MB) but odd for a zero phase 
difference (linear polarisation of the MB) as depicted in fig. |[ Thus, due to the parity properties of the 
vibrational states, linear polarization of the MB induces the maximum coupling between neighboring 
vibrational level whereas a phase difference of n/2 between the polarization components induces a 
coupling only between vibrational levels of the same parity. In the experiment, the relative phase 
between the two polarization components of the MB is controlled by an adjustable retardation wave 
plate. For a linear polarization of the MB making an angle a with respect to Ox (see fig. [l]), the 
coupling between \m = 3, n) — + \m = 2, n — 1) is 

V = ^i]U A Y ag (-^ - -^j sin (a) sin(6» Y AG)v / ^, = V R ^/n, (1) 

where Ai (resp. A2) is the detuning of the YAG beams with respect to the Cs D\ (resp. Dq) line, 
Ayag = Ai/3 + 2A2/3 and Uq = 4(Eodo) 2 / Ayag is the depth of the trap with linear polarizations. 
In this formula, only the component of order 1 in 77 has been retained, which is a good approximation 
in our situation where 77 = 0.16. For an angle a = 20°, the Rabi frequency of the coupling \m = 3, n = 
1) — ► \m = 2, n = 0) is 

Vr 

Q R = 2— = 2vr x 6 kHz. (2) 

h 

In order to measure the coupling due to the YAG beams, we first polarize the atoms in \F = 
3, m = 3) by an intense repumping pulse. We then measure the time evolution of the \F = 3, m = 3) 
population using a Raman transition to \F = 4, m = 4) as in fig. || (see also ref. (2^]). Figure |] gives 
the evolution of the population in \F = 3, m = 3). The magnetic field for this experiment is chosen 
so that the levels \m — 3, n) and \m = 2, n — 1) are degenerate and the angle of polarization of MB is 
a = 20°. We observe a damped oscillation. A damping with a time constant r ~ \/Tilo osc / (k^T)Vlii 
is expected as the coupling |m = 3, n) — > |m = 2, n— 1) is proportional to \fn. For a temperature of 
26 /iK, we expect r ~ 0.4^^. The Raman coupling, averaged over the population of the vibrational 
states, is obtained by fitting the first 40 [is of the oscillation with a sine function of amplitude 1. We 
obtain a Rabi frequency of 5 kHz in reasonable agrement with eq. |2| In this analysis we have supposed 
that all atoms contribute to the oscillations. However, this is not strictly the case. First atoms in 
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FIG. 4: Experimental determination of the YAG coupling \F — 3, m = 3, n) — > |_F = 3, m = 2, n — 1). 
The atoms are previously polarized in \F — 3,m = 3} by a strong pulse of the repumping beam. Crosses 
: evolution of the number of atoms transfered from \F = 3.m = 3) to \F — 4, m = 4) as a function of the 
delay time of the Raman probe pulse. Dashed line : fit of the first 40 /j,s with a sine function of amplitude 1. 



Squares : number of atoms transfered to 
\F = 3, m = 2} -> \F = 4, to = 3) (see fig. 



F = 4) on the transitions \F = 3, m = 3) — » \F = 4, m = 2) and 
3 and ref. P3J). We deduce a Rabi frequency ~ 5 kHz. 



\m = 3, n = 0) are not transfered to \m = 2). Second, atoms far from the center of the trap have a 
vertical oscillation frequency which differs from the central one. When the difference exceeds ~5 kHz, 
the coupling is no longer perfectly resonant. Both of these effects would lead to a Rabi frequency 
slightly higher than the estimation of 5 kHz. 



IV. COOLING OF THE VERTICAL MOTION 



A. Measured temperature 



The smallest vertical temperature we achieved with this sideband cooling method is about 0.56?iu;osc, 
which corresponds to about 83(5)% of the atoms accumulated in the ground state of motion. This 
temperature is measured either by time of flight absorption imaging or by using a Raman transition 
from \F = 3,m = 3) to \F = 4, m = 4) with counterpropagating beams. The vibrational sidebands 
are resolved (fig.[|) and from the height of the red sideband we deduce the population of the n > 
vibrational states Both measurements are in agreement at the 5% level. 

The measured 1/e cooling time constant is 2.5 ms. Such a cooling is obtained for an intensity of 
the repumping beam of 2.4 x 10 _2 / sa t, where 7 sat = l.lmW/cm 2 is the saturation intensity. This 
beam intensity corresponds to a calculated energy width of \m = 2) hV = h x 4.8 kHz and a linear 
polarization of the MB with a — 20° corresponding to a Raman coupling |m = 3, n = 1) — ► |m = 
2, n = 0) of h x 2.5 kHz. This ground state population of 83% is slightly less than the value (92(5)%) 
we previously obtained with sideband Raman cooling of unpolarized atoms involving the two hyperfine 
states \F = 3) and |F = 4) Jl|. 



B. Discussion 



The expected equilibrium temperature can be easily calculated in the Lamb-Dicke regime where 
77 1 if the equilibrium is set by non-resonant excitation from \m = 3, n = 0) to \m = 2), pj , p5| , p6| . 
The following simple argument gives the equilibrium temperature in the limit where the Raman 
coupling is much smaller than the lifetime r' of \m = 2). For 77 -C 1, we can neglect any change of 
vibrational level during the repumping process. Furthermore, the Raman transfer inducing a change 
of the vibrational level by k > 2, of order 77 fe , can be neglected compared to the Raman coupling with 
a change of vibrational level by 1. Thus, the two processes in competition that define the equilibrium 
are the process \m = 3,n = 0) — > \m = 2, n = 1) \m = 3,n = 1) and the resonant process 
\m = 3, n = 1) — > \m = 2, n = 0) => |m = 3, n = 0), where => indicates optical pumping. Note that 
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FIG. 5: Number of atoms transferee! to the hyperfirie ground state \F = 4) by a two-photon Raman pulse as 
a function of the detuning of the two photon transition from 9.192631770 GHz. The two Raman beams are 
counterpropagating in order to enable change of vibrational level. The magnetic fiel has been turned off before 
the Raman pulse so that all the transitions \F = 3,m) — > \F = 4, m') are degenerate. Solid line : spectrum 
taken after 20 ms cooling. Dashed line : spectrum taken before cooling. The labels indicate the change of 
vibrational level during the transfer. On the spectrum taken after cooling, contrary to the spectrum taken 
before cooling, the sideband corresponding to the transfer with a decrease of the vibrational level by 1 has 
deasapeared. This is expected if all the atoms are in the motional ground state \n = 0). 



the Raman transfer |m = 3, n = 0) — ► \m — 2,n — 0) is not allowed with linear polarization of the 
MB. In the limit where the Raman coupling is much smaller than the energy width T' of |m = 2), 
equaling the rates of these two processes gives a relative population of \m = 3, n = 1) of (r'/4tj osc ) 2 . 
The expected cooling rate is r coo i = fl^/T' . With our experimental parameters, this leads to a 
temperature of about 0.12?kj osc (population of |m = 3, n = 1) of 0.0002), a factor of 4.6 smaller than 
the obtained temperature. The predicted cooling rate is 3 x 10 4 s _1 (1/e cooling time of 30 fis). 

We discuss now some possible sources of discrepancy between this simple model and the experiment, 
although we find that they are not sufficient to explain the observed equilibrium temperature. First, 
the resonance of the transition \m = 3, n = 1) — > |m = 2, n — 0) is not fulfilled for all the atoms. 
Indeed, the vertical oscillation frequency of the atoms, which depends on their position in the trap, 
spreads over about 15% of the central 80 kHz frequency. If the transition is detuned by 6 > T', the 
expected steady-state population of |m = 3, n = 1) is (<5/(2w osc )) 2 . For 5 = 0.15w OS c, this population 
is only 0.005, still smaller than the measured one. However, the cooling rate is significantly modified 
by this effect; for instance, it is reduced by a factor 36 for a detuning of 0.15huj oac = 12 kHz. Thus, the 
spread of oscillation frequencies is likely to be the cause of the long cooling time observed. Second, 
non resonant Raman transfer is not the only process causing a departure from \m — 3, n = 0). 
Indeed, a component of cr_ polarization of the repumping beam would excite atoms in |m = 3). The 
measurement of the steady-state polarization without Raman coupling shows that the excitation rate 
of \m = 3) is smaller than O.Oir'. The heating produced by this excitation is smaller than or on the 
order of 0.01 x _E roc r'. For atoms resonant with the |m = 3, n = 1) — > \m = 2, n = 0), this would lead 
to a steady state population in the excited vibrational states of 0.003 much smaller than the measured 
one. However, atoms detuned by 12 kHz (the typical width of the distribution of oscillation frequency 
in the cloud) experience a smaller cooling rate and therefore the steady-state population of excited 
states increases to about 1% for these atoms. 

Finally, heating due to multiple photon scattering within the cloud is likely to bring a negligible 
contribution to the observed mean vibrational number because the optical density is only 2. 



C. Effect of the polarization of the MB 

The Raman transfer |m = 3, n) — > \m = 2,n— 1) is resonant when the Zeeman splitting is equal 
to the vibrational energy. Figure [(] (a) shows this resonant behavior of the cooling as a function of 
magnetic field. For a Zeeman splitting equal to 2lo osc , a transfer \m — 3, n) — > \m — 2, n — 2) is 
resonant. But with a linear polarization of the MB, such a transfer is forbidden (no coupling between 
levels of same parity as shown in fig. || of section [II C| ). In contrast, when we choose an elliptical 
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polarization of the MB, cooling on the second sideband is observed (fig. || (b)). This cooling is weak 
because the Raman coupling is very small for a small Lamb Dickc parameter. In this case, as expected, 
almost no cooling is observed on the first sideband. 
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FIG. 6: rms width of the vertical velocity distribution obtained after 30 ms of cooling as a function of the 
Zeeman splitting between adjacent Zeeman levels, vo = y^hmuj osc /2 is the rms width of the vibrationnal 
ground state. The polarization of the MB is linear in (a) and elliptical in (b). 



V. COOLING THE HORIZONTAL MOTION 



To cool the horizontal degrees of freedom with the vertical sideband cooling, a transfer of energy 
from horizontal to vertical motion is needed. This transfer can be provided by collisions between 
atoms and efficient cooling of the three degrees of freedom has been demonstrated in Q|. In our 
experiments, after 1 s of vertical sideband cooling, the temperature drops from Th = T v = 13 /iK to 
Th = 2.7 /iK and T v = 1.6 /iK as presented in fig. 0. For these data, the cooling rate of the vertical 
motion is much higher than the horizontal one. The intensity of the repumping beam is decreased to 
about 2.5 x 10~ 3 7 sa t, which is about 10 times smaller than the power that leads to the most effective 
cooling of the vertical motion. This is due to the heating of the horizontal motion experienced by 
the atoms even after reaching the steady state of the vertical motion. Indeed, the steady state of the 
vertical cooling is an equilibrium between heating processes like excitation by the er_ component of 
the repumper and the cooling itself. Spontaneous photons are emitted at a rate proportional to the 
intensity of the repumping laser and the spontaneous photons emitted in the horizontal directions are 
responsible for a heating of the horizontal motion. To achieve a heating smaller than the cooling rate 
due to the collisions, the power of the repumper should thus be small enough. Figure [?] also indicates 
that the cooling in the vertical direction is affected by the energy in the horizontal planes since the 
rms velocity in the vertical direction tends slowly to the ground state rms velocity as the horizontal 
temperature decreases. 

In the horizontal plane, we could achieve a temperature T = 2.8 /iK (ksT — 0.7tko osc ) which cor- 
responds to a mean kinetic energy per degree of freedom of 1.37kJ sc/4. For such a low horizontal 
temperature, the energy transfer between the horizontal and ver tical degree of freedom due to the 
collisions is very inefficient. Indeed, we explain in section VIII E , that for small horizontal tempera- 
tures, the rate of energy transfer due to collisions goes down exponentially with a factor e - 2h ^c/(k B T) 
as the temperature decreases. This factor, which is 0.5 for the typical initial horizontal temperature 
(/cbT = 3?kj osc ), is only 0.06 for the final temperature. Thus, even if the vertical equilibrium temper- 
ature of the sideband cooling T eq is much smaller than hu osc , we expect that the rate of the horizontal 
cooling decreases exponentially as ksT < tuv osc . Small external sources of heating would then cause 
the horizontal cooling to stop before T cq has been reached. This inhibition of the coupling between 
the vertical and horizontal degrees of freedom for temperatures smaller than hu osc is an important 
limitation of 3D cooling using ID sideband cooling and collisions. 

The highest phase space density we reached is hAq B = 1.3 x 10~ 3 , where Adb = hy / 2n/mkBT is 
the de Broglie wave length and n is the peak density. This is obtained after cooling to a temperature 
fc^Tft = k&T v — 4.3 /iK with about 450 atoms per micro-traps at a peak density n = 4 x 10 12 
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FIG. 7: Evolution of the rms velocity width in the vertical (+) and the horizontal (x) directions as a function 
of the cooling time, vo is the width of the ground state of motion in the vertical direction. Note, at very short 
cooling time the initial increase of the velocity spread in the horizontal direction. It is due to the spontaneously 
emitted photons required for the sideband cooling of the vertical motion. 



atoms/cm 3 . 



VI. A METHOD TO COOL FURTHER : CHANGE OF POLARIZATION ANGLE a 

To overcome the decrease of the cooling efficiency of the horizontal motion as the temperature 
decreases, we changed the parameters of the trap during the sideband cooling. More precisely, the 
ratio of the horizontal to the vertical oscillation frequencies is increased. Thus, if the change is done 
adiabatically with respect to the oscillation of the atoms, the ratio of the energy of the horizontal 
motion to the energy of the vertical motion will increase. More pairs of colliding atoms will then 
have an horizontal energy sufficient to populate the excited vibrational states after the collision. This 
change is realized by increasing the angle a of the polarization of the MB with respect to the horizontal 
direction, using a retardable wave plate on a time scale of about 10 ms. With a large angle a, the 
contrast of the interferences between the two trapping beams is reduced leading to a decrease of the 
vertical oscillation frequency by a factor ^/cos(a). Because the depth of the trap is also reduced, the 
horizontal oscillation frequencies are also reduced but by a smaller factor of + cos(a))/2. 

Experimentally, we decrease the vertical oscillation frequency from u> Zl = 80(3) to uj Z2 — 53(2) kHz, 
(decrease by a factor 0.67), corresponding to a going from 29° to 63°. The horizontal oscillation 
frequencies decrease only by a factor oj X2 /u) Xi = 0.88. At 53 kHz, sideband cooling still works well, 
after the corresponding change of the magnitude of the magnetic field. We also found that a decrease of 
the repumping power by a factor 3.3 was required. Figure || gives the time evolution of the horizontal 
and vertical velocity distributions before and after the change of the vertical oscillation frequency. 
More cooling after the change of the oscillation frequency is clearly visible. 

If we assume that the cooling for a large value of a produces ksTh = f3Tiuj Zl when returning 
the polarization to its original value this final horizontal temperature is given by fc^T/j/ '(Jiu) Zl ) = 
(3(uo Xl /uo X2 )(ui Z2 /uo Zl ). This way, an horizontal temperature that could not be obtained by cooling 
only with oj Zi because of the decoupling of the horizontal and vertical motion at low temperature 
could be produced. 
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FIG. 8: Two step cooling with a decrease of the vertical oscillation frequency after 500 ms. squares : rms 
velocity width along z. circles : rms velocity width along x. The dashed straight lines give the rms velocity 
width of the ground state of motion in the vertical direction. 
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FIG. 9: Width of the rms velocity in the horizontal (dashed line) and vertical (solid line) direction after 60 ms 
(a) and 500 ms (b) cooling time as a function of the number of atoms measured at the end of the cooling. 



VII. DEPENDENCE OF THE TEMPERATURE ON THE NUMBER OF ATOMS AND 

ATOM LOSSES 

Dependence of the measured temperature with the number of atoms has been observed as shown 
in fig. |^. At large atom numbers the temperature is higher. This can be due to heating by the 
reabsorption of spontaneous photons or by exothermic light-assisted collisions. 

We also found that this sideband cooling method is accompanied by losses of atoms : about 50% 
of the atoms have been lost in the first 500 ms of the cooling presented in fig. (?]. The corresponding 
loss rate of about 1.0(1) s _1 is significantly larger than the measured loss rate due to collisions with 
the back-ground gas (0.4(1) s _1 ). 

We believe that these losses are due to exothermic collisions assisted by the repumper light [p7| . 
However, a dependence of the loss rate on the density of atoms characteristic of such a 2-body loss 
mechanism has not been demonstrated. Indeed, due to the dependence of the temperature on the 
number of atoms, only small changes in the density have been achieved (~ 15%) and no conclusive 
results could be deduced. 

For the depth of our trap (about 130 /iK), the dominant inelastic light-assisted process is radiative 
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escape j28). We give below an estimate of the loss rate due to this process in our experiments, based 
on earlier calculations made by Julienne and Vigue |29|| . The loss rate can be written Kn ex , where 
n ex is the density of excited atoms. For a temperature of 100 fiK, a trap depth of IK, and a detu ning 
of the laser of — r, K ~ 0.2 x 10 -11 |§]. Because K scales with the depth U of the trap as U 5/6 [§9j, 
we expect K to be bigger by a factor 1.7 x 10 3 in our trap of depth 130 fiK. In our experiment, the 
rcpumping beam is on resonance and we could thus expect a change of the factor K . However, we 
believe that it will be of the same order of magnitude. The density of excited atoms after the first 
10 ms of strong vertical cooling, is mainly determined by excitation of the atoms in \m = 3) due to the 
er_ component of the repumper. With the experimental parameters and taking a component of bad 
polarisation of the repumping beam of 5%, n ex is about n x 0.6 x 10~ 4 , where n is the total density 
of atoms. With a mean density n = 5 x 10 11 at/cm 3 , the estimated loss rate is then 0.15 s -1 . This 
loss rate is of the same order of magnitude as the measured one. 

We note that, contrary to our findings, no atom losses were observed in ref|jj], although the density 
achieved in Q was 3 times higher than ours. 

VIII. COLLISIONAL PROPERTIES 
A. Introduction 

In the ground state of motion in the vertical direction, atoms are confined in a gaussian distribution 
with rms size Iq = 20 nm. On the other hand, the scattering length of cesium atoms which characterizes 
the collisional properties at low energy is, for the \F = 3, m = 3) state, negative with an absolute value 
larger than 60 nm j|, |l9| , |2 ^ | (The most recent calculations of ref||, based on the analysis of several 
Feshbach resonances of j22j give a = — 138 nm). An interesting question is then : are the collisional 
properties altered by the strong ID confinement as compared to the free space case ? With such a 
high absolute value of the scattering length, for the temperatures in our experiments, the collisional 
cross section for free atoms reaches its maximum value of 8ir/k 2 . This resonance behavior has been 
observed in weak traps by measuring thermalization times as a function of temperature ]l9| , |2p| , [30| . 
The quantity l/(rw Ttherm), where n is the density of atoms, which is expected to be proportional to 
the typical collision cross section a, was found proportional to X/T. This was interpreted as a zero 
energy resonance and a was found equal to its maximum allowed value of 8ir/k 2 . We give below an 
analytical derivation, in the classical limit (fc^T 3> tko osc ), of the thermalization time T t herm- 

In our trap, at temperatures much larger than the vibrational energy huj OSCl we expect the collisions 
to be 3D collisions. Indeed, in this case, the motion of the atoms can be understood as oscillations of 
wave packets of quite well-defined wave vector k. The collision of two wave packets occurs on a spatial 
scale l/k (square root of the cross section). For the typical k = ^/rnksT JK and for ksT ^> Tilo OSCi this 
spatial scale is much smaller than the spatial amplitude of the oscillation hk / (muj osc ) . 

However, at temperatures kgT < htu osc , D. Petrov and G. Schlyapnikov predict drastic changes of 
the collisional and thermalization properties, such as a change of the sign of the scattering length^] 
and an exponentially vanishing rate of the thermalization between the strongly confined motion and 
the transverse onesQ. As a test of these possible changes we have investigated the zero-energy 
resonance and searched for confinement-induced modifications of the thermalization rates at various 
temperatures. 

B. Measurements 

A situation out of thermal equilibrium is easily produced by the vertical sideband cooling method 
described above. As the cooling rate of the vertical motion is much larger than the collision rate, 
most of the atoms are quickly cooled to the ground state of the vertical motion, while the horizontal 
motion is cooled more slowly. At any cooling time, the horizontal temperature is higher than the 
vertical one. If the vertical cooling is stopped, the cloud of atoms relaxes towards equilibrium. We 
then measure the time evolution of the widths of the velocity distributions along the vertical direction 
and the horizontal direction. We checked that this relaxation is indeed due to collisions by measuring 
thermalization times for different numbers of atoms : Ttherm is found to be inversely-proportional to 
the number of trapped atoms. 
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The initial horizontal temperature is varied by changing the time duration of the sideband cooling 
phase before the measure of thermalization. Wc could thus measure thermalization times for initial 
horizontal temperatures varying from 20 fj,K to 4/iK. The maximum temperature is limited by the 
depth of the trap. Figure ^ presents many measurements of thermalizations corresponding to various 
initial temperatures. As the absorption imaging is destructive, each point in fig. ^ corresponds to a 
different experimental cycle lasting about 2 s. Fluctuations in the number of atoms for a given cooling 
time are about 2%. We check that 2E C:c + E Cz , where E Cm and E Cz are the mean kinetic energies along 
the vertical and horizontal axis of the camera, is conserved to within about 5% (see fig. |10|). This 
behavior is expected if the thermalization time is much larger than the oscillation time. Indeed, in 
this case, we expect that the energy is always equally distributed between the potential and kinetic 
energy and that the horizontal velocity distribution stays isotropic. 



C. Results 

Figure [OJ presents, as a function of temperature, the quantity l/(ri3£>tv ms Ttherm) where v rms is 
the initial (i.e. at the beginning of the thermalization) rms width of the horizontal motion and T t h e rm 
is the 1/e thermalization time deduced from an exponential fit to the evolution of the vertical velocity 
width. ri3D is the initial mean density seen by an atom. It is averaged over the different vertical 
micro-planes and the population of the different micro-planes is deduced from the measured gaussian 
vertical distribution of atoms and the lattice period. Furthermore we assume that the temperature 
is identical in all micro-planes, n^jj is then deduced from the temperature measurement and from 
the vertical and horizontal oscillation frequencies. The horizontal oscillation frequencies (of the order 
of 175 and 140 Hz) are measured through the heating of the atomic cloud produced by parametric 
drive |3l| with a precision of about 5%. 

The straight line in fig.[ll] corresponds to the expected behavior when kgT hto osc . It results 
from the 3D classical calculation presented below, where the scattering resonance is assumed but 
where the quantization of the vertical motion is not taken into account. Within our experimental 
accuracy and for the explored 4 /iK- 20 /iK temperature range, we have seen no modification of the 
scattering resonance induced by the strong ID confinement: our measurements are compatible with 
the classical calculation. They are also compatible with the full quantum calculation performed by 
Petrov et al. Q (dashed line in fig.pd]), which departs significantly from the classical calculation only for 
temperatures much below 4 /iK. The quantum calculation assumes an oscillation frequency of 80 kHz 
and a scattering length a of — 60 nm. As this scattering length is already larger than the size of the 
ground state wave function, the calculation does not differ much from the a = oo case. It also shows 
that for fcgT -C fiuj osc the thermalization time is expected to increase exponentially as 

_ 9m 1 fi"°.c 

Ttherm = TTT" e k B T (3) 

64ft n 2 D 

As we will show below, this behavior can also be understood assuming that the scattering resonance 
is still valid but taking into account the discrete structure of the vibrational energy in the vertical 
direction and parity conservation during the collision. 



D. Classical calculation of Tthcn 



We propose here a calculation of the thermalization time in the limit ks T 3> tiu; osc . In this condition, 
the motion of the atoms can be treated classically. The collisions are the same as for free atoms and 
can be considered as "point-like" because the extension of the cloud is much bigger than the square 
root of the typical cross-section. Finally we assume that a = 8ir/k 2 . 

We are interested in the description of the thermalization between the horizontal motion and the 
vertical one. Thus we consider the evolution of 

AE = E Z ~ X -E p , (4) 

where E z and E p are the kinetic energies of the vertical and horizontal motions. AE vanishes at 
thermal equilibrium. If we neglect the anharmonicity of the trapping potential, only collisions can 
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FIG. 10: Free evolution of the rms width of the velocity distribution along the vertical direction (stars,v z ) and 
along the horizontal direction x (open square, v x ), after switching off the sideband cooling for different initial 
conditions. The solid line is an exponential fit of the evolution of v z and the 1/e time is noted on the graphs. 
On the second and fourth graph, the evolution predicted by the calculations of D. Petrov and G. Shlyapnikov for 
the same initial velocity distribution and for the same initial density are shown in dashed lines j^/. The second 
column shows y (vl + 2v%)/3, the velocity corresponding to the mean kinetic energy per degree of freedom. 
The third column gives the number of atoms. Each point is the average of 20 successive measurements. 



change AE. A collision where two atoms of velocities Vi and V2 end up with the respective velocities 
V3 and V4 changes AE 1 by 

^12-34 = \m (v 2 Z3 + v 2 Z4 - v 2 Zl - v 2 Z2 - i (v 2 p3 + v 2 p4 - v 2 pi -v 2 2 )Y (5) 
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FIG. 11: Comparaison between experimental results of thermalization and theoretical calculations in the tem- 
perature interval [4 [iK, 20 p,K\ . The solid line corresponds to the 3- dimensional behavior (eq. The dashed 
line results from the quantum calculations of D. Petrov et al. (see text). Two typical error bars at high and 
low temperatures are shown. . 



where the indices p and z refer respectively to the horizontal and vertical component of the velocity. 
Energy conservation implies 

c ^^/2 2 2 2 \ / r>\ 

<>12^34 = 2 2™ + V Z4 ~ V zi ~ V z 2 )- ( 6 ) 

If we note /(r, v) the phase space distribution, the rate of change of AE is 
dAE 



dt 



J d 3 rj J d 3 v 1 d 3 v 2 d 3 v 3 d 3 v 4 <5i2- ( 34Wi2^ 3 4 (/(r, vi)/(r, v 2 ) - /(r, v 3 )/(r, v 4 )) , (7) 



where Wi2^34 = W^^xz are the collision rates. The factor 1/4 cancels the fact that each collision 
process (1 <-> 3, 2 <-> 4) is counted 4 times. 

Changing variables in the integral, using the center of mass velocity and the angle of deviation in 
the center of mass frame, W is 

^12^34 = |v 2 - vi|-^. (8) 
An 

with (7 = 87i7fc 2 , where k — |v 2 — Vi|m/(2?i) is the wave vector of the relative motion. 

Equation [7] alone does not give the evolution of the distribution /. In general, one has to solve the 
Boltzmann equation. But, in the collisionless regime where the time between two collisions is much 
longer than the oscillation period of the atoms in the trap, the problem can be simplified. In this 
regime, we do not expect oscillations of the cloud and we can assume that, for both the vertical and 
horizontal motion, the distribution / corresponds to a Boltzmann distribution with temperatures T z 
and T p which may be different. 

With such an ansatz for /, and noting that 

AE = Nk B (T z ~T p ), (9) 

eq. [7] gives a differential equation for T z — T p . For small initial deviations from equilibrium, we can 
linearize the equation in T z — T p . With this approximation, the evolution is exponential with a 1/e 
time constant 

_ 150F 1 

Jthcrm - 7, 1 1 ( iU ) 

k BT~ \l m 

where Tq is the equilibrium temperature and n is the mean density. 
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We thus find that 



-cr(v rms ) = — — — , (11) 



nv rms T theim 150F rms 15 mk B T ' 



where a(v rms ) is the cross section for a relative velocity v rms — ^/kj^To/m. 

As expected, the quantity •=- — L is inversely proportional to the temperature. We have thus 

obtained the numerical factors entering in eq. [ll], which is plotted in fjg.|ll|. 

E. Inhibition of thermalization when ksT <C hui osc 

We now give a simple physical interpretation of the exponential decrease of the thermalization 
rate between the horizontal and vertical degrees of freedom previously found in gj. We show that 
this behavior is due to the quantization of the energy levels in the vertical direction and to parity 
conservation during the collision process. 

The collision rate is 



coll 



(12) 



where n is the three dimensional density, v is the typical velocity of the atoms and a is the typical 
cross section. For fcgT <C Tiu osc , the horizontal velocity of the atoms is much smaller than the typical 
vertical velocity yjhuo osc / 1 (2m) and the 3D density is n — niD^JImj (uj osc h), where n2D is the 2D 
density. Thus the collisional rate satisfies 

r co ,i^^. (is) 

m 

■^coii * s an estimate of the thermalization time between the two horizontal degrees of freedom, in 
agreement with the rigorous calculation of ref. || as long as the temperature is not extremely small 

ID- 

On the other hand, the thermalization between the horizontal and the vertical degrees of freedom 
actually takes a much longer time. Indeed, the energy levels are quantized in the vertical direction. 
In order to transfer energy to the vertical motion during a collision, the horizontal kinetic energy of 
the atoms has to be at least hio osc before the collision. In fact, because of parity conservation, only 
the vibrational state \n = 2) of the relative motion can be populated after the collision M. Thus, 



the horizontal energy of the pair of colliding atoms has to be at least 2Uuj osc . Figure 12 presents 
the probability distribution n(E) of the kinetic energy of the relative horizontal motion of a pair 
of atoms, for a temperature of 3?ia;osc and a temperature of 0.5Hlo osc close to the value achieved in 
our experiments. Because the density of states is constant in 2 dimensions, at a temperature T, 
n(E) — 1 / ' (kBT)exp(-E / '(fc^T)) . Thus, the probability that a colliding pair has an energy large 
enough to populate the excited states of the vertical motion is 

p = e "b t . (14) 

Among the collisions transferring energy to the vertical motion, we can neglect the exponentially 
small amount that will transfer an energy larger than 2huj osc . Thus, the rate of change of the vertical 
energy is 

— - ~ 2huj osc ^^e "bt . 15 
at m 

This exponential factor is responsible for the drastic inhibition of the cooling of the horizontal motion 
near k^,T ~ Tiuj osc as already discussed in section 

On the other hand, the increase of the vertical energy needed to reach the thermal equilibrium is, 
neglecting the change of horizontal temperature, 

AE Z = ftw osc e"^. (16) 



Thus, the thermalization time Xthcrm is of the order of 



1G 



1 dE z n 2D h 

e k B T (17) 



Ttherm eft 

This formula is identical, within a numerical factor, to the formula of ref.H using a quantum 
treatment of collisions in 2D. It is interesting to note that both dE z /dt and AE Z are exponentially 
small but it is the difference in the exponents that explains the exponential factor of T t h crm . 




FIG. 12: Probability distribution of the relative kinetic energy of the horizontal motion of a pair of atoms for 
a temperature of 3huj osc / ks ( initial conditions in our experiments) and a temperature of 0.5hu! osc /kB ( end 
of ID sideband cooling. 



The most efficient way to detect this exponential dependence of Tiherm would be to compare the 
thermalization time between the two horizontal degrees of freedom (~ Tiri2D / fn) to the thermalization 
time between them and the vertical motion. However when the initial horizontal temperature is 
fceT < fr^osc, even if the initial vertical temperature is ~ 0, the relaxation towards equilibrium 
corresponds to an exponentially small change of the horizontal temperature and an exponentially small 
population in the excited vibrational states. The observation of these small changes would require a 
very sensitive measurement of the temperature or of the population of the vertical vibrational states. 
Therefore the easiest choice of parameters to detect this 2D behavior is to operate near k^T ~ Huj osc . 



IX. SUMMARY 



In this paper we have investigated ID sideband cooling of cesium atoms in a far-detuned optical 
lattice. This cooling is particularly efficient because, for each spontaneously emitted photon, an energy 
of order huj osc is removed from the system. Atoms are spin-polarized in F=3 and, at high density, 
elastic collisions enable 3D cooling of the sample through ID sideband cooling. We have produced 
atomic samples in which more than 80% of the atoms are in the vibrational ground state of motion 
in ID and with a transverse temperature of 0.7 hw osc /k-Q = 2.8 /iK. This realizes a quasi-2D cold gas. 
Limitations of this cooling method have been identified. First, light-induced atom losses have been 
observed. Second, cooling of the weakly confined degrees of freedom by collisions looses its efficiency 
around k^T ~ Tiuj osc . We explained this effect by an exponential slowing down of the energy transfer 
between the weakly and strongly confined degrees of freedom. A dynamic method involving changes 
of the high oscillation frequency via polarization rotation of one of the trapping beams has been 
proposed and implemented. We have shown that the zero energy scattering resonance of cesium was 
essentially unaffected in the temperature range 4/iK— 20/xK. Within our experimental accuracy, this 
is in agreement with theoretical calculations. Further studies could concentrate on collisions affecting 
only the motion in the weak confinement directions where dramatic modifications of collisions occur 
at k^T < 0.1?iw osc [j2|, |6|. Promising systems for these studies are Bose- Einstein condensates loaded in 
a far-detuned trap associated with the possibility of manipulating the scattering length via Feshbach 
resonances. 
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